Astronomy & Astrophysics manuscript no. 
(will be inserted by hand later) 



Multi-wavelength spectrophotometry of EX Hydrae * 

S. Eisenbart, K. Beuermann, K. Reinsch, and B.T. Gansicke 

Universitats-Sternwarte, Geismarlandstr. 11, D-37083 Gottingen, Germany 
Received October 1, 2001 / Accepted November 23, 2001 

Abstract. We present phase-resolved infrared and optical spectrophotometry of the intermediate polar EXHya 
supplemented by archival ultraviolet data. The spin-modulated emission from the accretion funnel and the emission 
from the accretion disk or ring contain substantial optically thin components. The white dwarf dominates the 
unmodulated flux in the ultraviolet and is identified by numerous absorption lines. Metal absorption in the 
accretion curtain may add to the observed spectral features. The secondary star is of spectral type M4±l and 
is detected by its ellipsoidal modulation. We derive a distance of 65 ± 11 pc which makes EX Hydrae one of 
the closest cataclysmic variables with a known distance. The luminosity derived from the integrated overall 
spectral energy distribution is 3 x 10 32 ergs -1 . The accretion rate of 3 x 10 15 gs _1 (for an 0.6 A^q white dwarf) 
is in reasonable agreement with the rates expected from angular momentum loss by gravitational radiation and 
from the observed spin-up of the white dwarf. 

Key words. Stars: individual: EX Hydrae - cataclysmic variables - intermediate polars - accretion 



1. Introduction 

EX Hya was discovered by Kraft (1962) and quickly 
recognised as an eclipsing system with an orbital pe- 
riod P rb = 98 min. A second prominent periodicity 
Pepin = 67 min was interpreted as the rotation period of 
the white dwarf (Vogt et al. 1980, Kruszewski et al. 1981) 
which led to the intermediate polar model of EX Hya 
(Warner 1983). X-ray emission in the quiescent state was 
discovered by Watson et al. (1978), a hard X-ray eclipse 
was first seen by Beuermann & Osborne (1985, 1988). 

In this paper, we report new phase-resolved infrared 
and optical spectrophotometry. We have identified the sec- 
ondary star in the infrared and use its -fT-band flux to de- 
termine the distance to EX Hya, a key parameter for the 
discussion of its luminosity and accretion rate. The derived 
distance is smaller than thought previously, requiring a 
current accretion rate of only 3 x 10 15 g s _1 (for an 0.6 M Q 
white dwarf), consistent with the rates expected from 
gravitational radiation and implied by the observed spin- 
up of the white dwarf (Hellier & Sproats 1992). Combining 
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our data with archival IUE and HST spectrophotometry, 
we determine the contributions to the overall flux dis- 
tribution from the rotating funnel, the accretion disk or 
ring, and the white dwarf. The short-wavelength ultravi- 
olet spectrum contains metal absorption lines which are 
a clear signature of the white dwarf and/or absorption 
in the accretion funnel which surrounds the white dwarf. 
Based on these findings, we discuss the overall energy bal- 
ance and present an internally consistent description of 
EX Hya as an intermediate polar. 



2. Observations and Data Analysis 

2.1. Observations and Archival Data 

We have observed EXHya in February 1997 during three 
subsequent nights in the near infrared (0.8-2.4 fim) with 
the Infrared Spectrometer (IRS) at the 4-m telescope of 
the Cerro Tololo Interamerican Observatory (CTIO) and 
during one nig ht in the optical (3600-10200A) with the 
ESO Faint Object Spectrograph and Camera 2 (EFOSC 2) 
at the ESO/MPI 2.2m telescope at La Silla. We com- 
plemented these data by archival time-resolved ultravi- 
olet spectrophotometry taken with HST and IUE in 1995. 
Table summarizes the relevant information. During all 
observations EXHya was in a quiescent state. For com- 
parison and calibration purposes in the near infrared and 
optical ranges, spectra of several M dwarfs with spectral 
types between dM2 and dM5.5 and of an AO standard star 
were taken. 
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Table 1. Summary of time- resolved observations used for the present analysis, tint and t rcs are the exposure times 
and and the effective time resolution. The values quoted for the HST FOS data refer to time-averaged spectra (see 
text). 



Telescope 


Instrument 


A 

(A) 


date / UT 


tint 

(sec) 


trcs 

(sec) 


Observer 


CTIO 4-m 


IRS 


8200-24200 


Feb 22/23 1997 05:55-09:55 
Feb 23/24 1997 07:35-09:20 
Feb 24/25 1997 07:00-09:55 


60 


65 


K. Reinsch 
S. Eisenbart 


ESO 2.2-m 


EFOSC2 


3600-10200 


Feb28/Marl 1997 06:50-08:50 


240 


- 360 


K. Reinsch 
S. Eisenbart 


IUE 


SWP 


1150-1980 


Jun 23 1995, 13:50- Jun 24 1995, 20:30 


600 


~ 2400 


K. Mukai 


HST 


FOS 


1160-1600 


Feb 12 1995 04:05-14:20 
Feb 15 1995 12:30-22:40 


3.34 


~ 13.3 


S. R. Rosen 



2.1.1. IRS 

The IRS at CTIO is a cooled grating spectrometer 
equipped with a two-dimensional InSb detector array 
which allows spectral resolutions of A/AA « 150. ..3000 
in the wavelength range 0.8 — 5 /zm (Depoy et al. 1990). 
When used by us, it was upgraded compared to Depoy et 
al. We chose a setup which allowed to observe the whole 
wavelength range of 8300 — 24200 A with a spectral reso- 
lution of A/AA ~ 560 and a time resolution of 65 s. The 
seeing was 1.0"-1.3" during the first and 1.2"-1.8" during 
the second and third night. A fraction of the observation 
time in nights 2 and 3 was affected by clouds and the spec- 
tra taken during these times were discarded. All spectra 
were obtained through a 2" aperture. 

A total of 222 useful spectra of EXHya were taken, 
the majority during the first night. In order to correct 
the spectra for water vapor absorption in the infrared, 
we observed the nearby G5 star SAO 181198 at regular 
intervals. The NIR spectrum of this star has virtually no 
intrinsic spectral features apart from the hydrogen lines 
and the hourly spectra taken of this star allow a correction 
for H2O absorption, at least on this time scale. 

The M dwarfs observed were G1382 (dM2), G1443 
(dM3), G1402 (dM4), G1285 (dM4.5e), and G1473AB 
(dM5.5e). As for EXHya, spectra of nearby G0-G5 stars 
were taken for correction purposes. 

2.1.2. EFOSC2 

The optical spectra were taken a few days after the in- 
frared observations when EX Hya was practically at the 
same brightness level. The sky was clear with a seeing of 
1.7". We observed EXHya through a 2" slit with low spec- 
tral resolution (EFOSC2 grating Gl, A/AA ~ 150) cov- 
ering the wide wavelength range from 3600 to 10200 A. A 
total of 20 spectra were taken with an integration time of 
4 min and an effective time resolution of 6 min. They cover 
1.2 orbital or 1.7 spin cycles and include two eclipses. 



2.1.3. IUE 

From the IUE data archive, we selected a continuous set of 
45 SWP spectra (A = 1150-1980 A) obtained by K. Mukai 
over 1.3 days in June 1995. They are the same spectra used 
by Mauche (1999) in his analysis of EX Hya. The spectra 
were taken through the large aperture with an exposure 
time of 600 s (A0 6 7 = 0.15). For further description of the 
data including a list of the phases of the individual spectra 
see Mauche (1999). The SWP spectra were supplemented 
by long-wavelength spectra which are effectively not spin- 
resolved, however. 

2.1.4. HST 

From the HST archive, we obtained two sets of spectra 
taken by S. R. Rosen with the Faint Object Spectrograph 
(FOS) on Feb 12 and Feb 15, 1995. The spectra were ob- 
tained with the G130H grating and the 4.3" aperture cov- 
ering the wavelength range 1160-1600A at a FWHM res- 
olution of 1.5 A. There are 14 data sets of 120 spectra each 
which extend over orbital phases (j)g 8 = —0.15 to +0.20 
and thereby cover 14 eclipses. The data were taken with 
the FOS operated in rapid readout mode at an integra- 
tion time of 3.34 s, resulting in an effective time resolution 
of 13.3 s. In order to improve the S/N ratio of the in- 
dividual spectra, we averaged 12 subsequent spectra and 
thus reduced each of the 14 data sets to 10 spectra with 
a S/N ~ 8 and an effective time resolution of 160 s. We 
note that the eclipse in EX Hya of about 3 min full width 
is seen in the original data, but is lost at the reduced time 
resolution. 

2.2. Data analysis 

The IUE and HST data required no further reduction, 
except for the removal of small spectral sections in the 
IUE data which were affected by reseau marks. The optical 
spectra were processed with standard MIDAS routines. 
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The data structure and the intense sky background of the 
infrared spectra required special routines which one of us 
(SE) developed and integrated into MIDAS. 

Absolute flux calibration of near IR spectra is diffi- 
cult because suitable spectrophotometric standard stars 
are rare. To overcome this problem, we took spectra of 
several A-type photometric standard stars. The near IR 
spectra of such stars are blackbody-like and almost fea- 
tureless except for H and He absorption lines. This al- 
lows us to construct the intrinsic spectrum of each star as 
a blackbody for the effective temperature of its spectral 
type and to reproduce observed photometric JHK mag- 
nitudes by appropriate least-squares scaling. In all cases, 
the colours of the constructed spectra agreed with the lit- 
erature values to better than ±0.05 mag. Finally, the in- 
dividual spectral response functions derived from the con- 
structed and the observed spectra agreed with each other 
to within a few percent over the whole wavelength range. 
Using this procedure, the mean atmospheric absorption in 
the object spectra is already accounted for in the spectral 
response function. The remainder was removed using the 
spectra of the 'monitor' stars mentioned in section |2.1.l|. 



3. Light curve analysis 

EXHya shows variations on the 67-min spin period and 
on the 98-min orbital period. For the light curve analy- 
sis, we used the orbital and spin ephemerides of Hellier 
& Sproats (1992), including the sinusoidal term in the 
former. Both ephemerides are given in Julian Days and 
refer to the barycentric dynamical time definition. Phase 
098 = indicates the center of the optical (and X-ray) 
eclipse and phase 067 = the maximum of the optical (and 
X-ray) flux during the spin cycle. We confirm the finding 
by other authors (Hurwitz 1997, Mauche 1999; Mukai et 
al. 1998) that the orbital ephemeris needs an update (the 
eclipse appears at 0g§ = —0.02 instead of at 0.00. This 
small deviation is not of relevance for the present study. 

The flux-calibrated time-resolved spectra allow an 
analysis of the spectral variations as a function of phase 
and a light curve analysis for finite bands in wavelength. 
We start here with the light curve analysis. 

From the individual spectra, we extract the time de- 
pendence of the total fluxes in broad bands (dominated 
by the continuum) and of the fluxes in the prominent 
emission lines above the continuum. The optical and IR 
broad bands approximately agree with Bessell UBVRI and 
Johnson I JHK, respectively, while the ultraviolet bands 
were, in part, chosen to be free of emission lines. We fit 
the individual data trains of the continuum and the emis- 
sion line fluxes with model light curves of the form 

F(t) = C + Bt + A 67 cos 2tt (0 67 (t) + A0 67 ) 
+A 98 cos 2tt (^gg (t) + A0 98 ) 
~A ig cos 2tt (0 49 (O + A M (1) 

where the term B t accounts for a slow variation of the 
mean flux over the observation time, Aq 7 , A98, and ^49 



are the amplitudes of the modulations at the spin, the or- 
bital, and half the orbital period, 067 , 98 , and 049 are 
the corresponding phases, and A0 6 7, A0 g8 , and A04g the 
phase shifts relative to the ephemerides of Hellier & Spoats 
(1992). (Note that a phase shift A04g = 0.02 corresponds 
to a shift by 0.01 orbital periods). The only significant 
modulation in the ultraviolet occurs at the spin period. 
Spin modulation is also dominant in the optical and our 
data train is too short to obtain meaningful parameters for 
modulation at the other periods. Only the longer time cov- 
erage in the infrared allows to deduce the modulation at 
all three periods provided for in Eq. 0. Tables | and | list 
the least-squares fit parameters. Of the emission lines, P7 
is covered, too, but can not sufficiently well be separated 
from HeiA 10830. The mean visual and if -band fluxes in 
Table || correspond to mean magnitudes of EX Hya of 
V = 13.43 and K = 11.77 at the time of our observations. 

3.1. Spin modulation 

Figure [l| (upper panels) show examples of the IUE fluxes 
fitted with line 1 of Eq. (Q): (i) the total flux in the 1250- 
1950A range and (ii) the emission line flux of CivA1550. 
Both fits benefit from including the term B t in Eq. (|I]) . 
We conclude that there is a gradual increase of the ultra- 
violet flux over the time of the observation on top of the 
spin modulation. The folded data, after subtraction of the 
linear term, A + Bt, are shown in the lower panels. Tables 
and I list the mean spectral fluxes /, the integrated line 
fluxes F, the corresponding amplitudes A e7 , phase shifts 
A067, and the corresponding fractional modulations, e.g., 
A^67 = 2 Ae 7 /f max , where / ma x is the mean spectral flux 
at spin maximum. Me7 — 1 corresponds to 100% mod- 
ulation and M.Q7 > 1 for SimA1300 implies that the line 
which is in emission at spin maximum turns into absorp- 
tion at spin minimum. Allowance for a time-dependent 
amplitude would further improve the fit in Fig. |l|, but com- 
plicate the definition of the modulation amplitude. The 
HST data yield very similar spin modulation but do not 
require the term Bt. The modulation of the line fluxes 
exceeds the modulation of the continuum in the HST and 
IUE data as well as in the optical data (Tables || and |^). 

The infrared data continue the decline in the amplitude 
of the spin modulation seen in the optical bands. Fig. 
shows the result for the JT-band (left top panel) and the 
P/3 line (right top panel). Optical depth effects in the in- 
frared emission line region are important as indicated by 
B7. This line displays practically zero spin modulation of 
the integrated line flux (Table ||) , but varies in width such 
that broad emission lines with deep absorption cores ap- 
pear in the spin-modulated spectral flux (spin maximum 
- spin minimum) . Fig. [| (bottom panel) shows B7 and the 
absorption cores of the higher Brackett lines up to B13. 
Such effects are still present in P/3 and decrease in the 
higher Paschcn lines. 
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Fig. 1. Fits of line 1 of Eq. (|l|) to sample IUE light curves of EXHya. Left panel: total flux in the 1250-1950 A band, 
right panel: CivA1550 emission line flux. Top panels: Original time series and fitted model, bottom panels: data with 
the linear term A + B t subtracted and folded over the 67-min spin phase. 
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Fig. 2. Left panels: Spin and orbital modulation of EXHya in the IR (quasi JT-band: 20000 - 24000 A). The data 
points around the eclipse (grey) were excluded from the fits. Right panels: Same for the Paschen P/3 line. The fit in 
the lower right panel represents the sum of the modulations at the orbital and half the orbital period. 



3.2. Orbital modulation 



In our optical spectrophotometry, the orbital eclipse oc- 
curs slightly early at 0g§ = 0.98. There is an indication 
for a shallow maximum near 0g 8 ~ 0.8 as shown to be 
present by Hellier et al. (2000, see also Siegel at al. 1989). 
The light curves of the Balmer line fluxes are flat and show 
a weak broad eclipse, as expected if part of the emission 
line fluxes originates from the accretion disk or ring. 



Our infrared continuum data (lower left panel of Fig. j2|) 
show the orbital eclipse superposed on substantial modu- 
lation at half the orbital period (data points covering the 
eclipse are displayed in gray). The infrared eclipse is more 
pronounced at spin maximum than at spin minimum, but 
is present at all spin phases. At the same time, the infrared 
eclipse is much wider than the eclipses at all other wave- 
lengths, indicating that a much more extended structure 
than just the accretion funnel is partially eclipsed. The full 
width of about 10 min corresponds to an eclipsed object of 
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Table 2. Mean spectral fluxes / in box-car shaped wave- 
length intervals and amplitudes of modulation at the spin 
Aqy in units of 10 -16 ergcm~ 2 s _1 A _1 . A(/>67 is the shift in 
the spin maximum relative to the cphcmcris of Hellier & 
Spoats (1992) and 7^67 the modulation amplitude defined 
in the text. The lower section gives the corresponding in- 
formation for the modulation at half the orbital period in 
the infrared. 



Table 3. Mean integrated line fluxes F and modulation 
amplitudes A§-j in units of 10 -16 ergcm _2 s _1 . A067 and 
Mq7 are as in Table |[ The lower sections give the cor- 
responding information for the modulation at the orbital 
and half the orbital period for the Cerro Tololo infrared 
data. 



Emission line F A§i A067 Me,i 



AA (A) / A 67 A067 Mm 



a) Hubble Space Telescope: 



1255- 


1285 


1480 


319 ±24 


-0.010 


0.355 


1350- 


1365 


1590 


276 ± 22 


±0.01 


0.296 


1350- 


1380 


1630 


304 ± 23 


±0.005 


0.314 


1425- 


1525 


1500 


221 ±19 


±0.005 


0.256 


1570- 


1600 


1360 


261 ±17 


±0.015 


0.323 



b) International Ultraviolet Explorer: 



1255- 


1285 


1450 


264 ± 40 


±0.05 


0.308 


1350- 


1380 


1560 


356 ± 60 


±0.015 


0.372 


1425- 


1525 


1440 


211 ±52 


±0.02 


0.256 


1575- 


1625 


1220 


250 ± 49 


±0.035 


0.341 


1680- 


1780 


1030 


262 ± 35 


±0.025 


0.407 


1875- 


1975 


1030 


283 ± 39 


±0.025 


0.430 



c) ESO EFOSC2: 



3600- 


4000 


539 


140 ± 18 


±0.03 


0.411 


3900- 


5000 


286 


59± 11 


±0.035 


0.342 


5000- 


6200 


155 


26 ±6 


±0.05 


0.287 


5800- 


8000 


121 


22 ±4 


±0.04 


0.312 


7500- 


10000 


75 


11.4±2.4 


±0.04 


0.265 



d) Cerro Tololo IRS: 

9650-10650 52.0 5.28 ±0.32 ±0.035 0.184 

11000-13500 37.6 3.31 ±0.22 ±0.03 0.162 

14500-18000 18.5 1.34 ±0.10 ±0.03 0.135 

20000-24000 8.3 0.47 ±0.05 ±0.025 0.107 



AA (A) / A 49 M49 



a) Hubble Space Telescope: 



cm 1175 


13000 




1 n ri/i c 


u.yo4 


Lya 


27300 


13000 ± 1000 


±0.035 


0.646 


NV 1240 


10300 


6990 ± 480 


±0.005 


0.807 


Si III 1300 


6190 


7720 ± 260 


±0.04 


1.11 


C II 1335 


9400 


4870 ± 390 


±0.035 


0.682 


Si IV 1400 


23300 


17300 ±660 


±0.02 


0.853 


C IV 1550 


56400 


37500 ± 1400 


±0.00 


0.798 


b) International Ultravi 


nlpf W/TTi] nrpT' 






Si IV 1400 


21300 


14100 ± 1800 


±0.015 


0.796 


C IV 1550 


51600 


39400 ± 3000 


±0.005 


0.865 


Hell 1640 


6740 


8490 ± 980 


±0.045 


1.11 


A1III 1860 


4960 


4560 ± 720 


±0.030 


0.958 


c) ESO EFOSG2: 








Hell 4686 


1413 


453 ± 58 


±0.025 


0.484 


H/3 


11288 


5135 ±350 


±0.01 


0.625 


He I 5875 


2288 


1114 ±73 


-0.005 


0.654 


Ha 


10113 


3463 ± 238 


±0.01 


0.510 


d) Cerro Tololo IRS: 








PS 


2060 


412 ±14 


±0.020 


0.333 


P/3 


2080 


334 ± 15 


±0.020 


0.277 


B7 


585 


8±5 


±0.000 


0.026 


Emission line 


F 


A 9S 


A098 


M 9 S 


P6 


2060 


139 ±17 


±0.330 


0.126 


P/3 


2080 


157±17 


±0.275 


0.140 


B7 


585 


40±8 


±0.240 


0.130 



9650-10650 52.0 2.39 ±0.34 -0.005 0.088 

11000-13500 37.6 2.16 ±0.23 -0.02 0.109 

14500-18000 18.5 1.21 ±0.11 -0.03 0.123 

20000-24000 8.3 0.56 ±0.05 -0.04 0.127 



about 3 x 10 10 cm diameter centered on the location of the 
white dwarf, quite likely the accretion disk or ring filling 
about 80% of the Roche lobe. A shallow wide eclipse was 
also seen in the optical data of Siegel et al. (1989) in ad- 
dition to the narrow eclipse of the accretion funnel. In our 
infrared data, the time resolution is not sufficient to dis- 
entangle the funnel and disk contributions to the eclipse. 
The central depth of the mean iT-band eclipse is about 
1.2 x 1CT 16 ergcm~ 2 s~ 1 A" 1 or 1.8 mjy. 

The infrared emission lines display a different orbital 
modulation and a weaker eclipse. Fig. || (lower right panel) 
shows P[3 as an example. The light curve is practically 
flat in the interval 4>q& — 0.5 — 1.0 and shows a hump at 



Emission line 



F 



A0 4 



M 4 



PS 

P/3 
B7 



2060 
2080 

585 



48 ±15 
100 ±15 
30 ±5 



±0.025 0.045 
±0.015 0.092 
-0.025 0.098 



098 = 0.25 which is phase shifted with respect to the blue 
continuum hump (Hellier et al. 2000) by roughly 180°. The 
fit shown along with the data includes the 98-min and the 
49-min variations which represent the first two terms of 
a Fourier expansion ^2 A n cos 2n n 4>g& of the orbital light 
curve. We interpret the maximum of the Paschen line 
flux as due to the inner illuminated side of the bulge, 
whereas the blue continuum originates from its outside 
heated by the interaction of the stream with the disk. This 
bulge is a lso detected by its ultraviolet line absorption (see 
Sect. below). 
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3.3. Ellipsoidal modulation of the secondary star 

In the infrared, the principal orbital modulation of the 
continuum is a double wave of 49 min period. There is 
no significant modulation at the orbital period. The mini- 
mum at <frgs = 0.98 coincides with the center of the eclipse 
to within 0.01 in <frgg. Fig. || (left lower panel) shows the 
variation of the X-band flux, with the other infrared bands 
behaving similarly. Table || lists the relevant fit parame- 
ters for our four infrared bands. The modulation depth is 
M. = 0.13 in K and decreases to 0.09 in the quasi-/ band. 

The double wave looks like the ellipsoidal modulation 
of the secondary star. While such a modulation is ex- 
pected, studies of other CVs show that other light sources 
can mimic the ellipsoidal modulation. In the dwarf nova IP 
Peg, part of the double hump structure likely occurs in the 
accretion disk (Froning et al. 1999) and in the magnetic 
system AR UMa, beamed cyclotron emission may con- 
tribute (Howell et al. 2001). Model calculations of Kube 
et al. (2000) for magnetic CVs demonstrate that a double 
hump can also result from the accretion stream. Typically, 
however, the energy distributions of the additional light 
sources differ from that of the secondary star and phase- 
resolved spectrophotometry should be able to disentangle 
their contributions. E.g., in low-state spectrophotometry 
of the magnetic system UZ For, the ellipsoidal modula- 
tion in the near infrared can safely bex separated from 
the bluer cyclotron source (Schwope et al. 1990). 

Considering these examples, we accept that some of 
the observed double hump modulation in EX Hya may 
originate from the accretion disk, but assume also that 
this contribution is less important than in other CVs as, 
e.g. IP Peg. The reason is that EX Hya displays no equiv- 
alent to the double hump in the optical continuum (Siegel 
et al. 1989, Hellier et al. 2000) nor in the optical and in- 
frared emission lines. Also, the accretion disk in EX Hya 
is not well developed (King & Wynn 1999). Another ef- 
fect which could affect the observed double humped light 
curve is the eclipse of the secondary star at superior con- 
junction (09g ~ 0.48) by the accretion disk. If the disk 
is largely optically thin, however, as suggested below, the 
magnitude of this obscuration would be small. In order 
to account for any other contribution to the 49-min in- 
frared modulation in EX Hya, we take the amplitude of 
the ellipsoidal modulation to be A c u = aA^g with a < 1, 
where A e ii is the amplitude of the ellipsoidal modulation 
at = 0.5. In estimating a, we refer to the analysis 
of the ellipsoidal modulation in IP Peg by Froning et al. 
(1999) which yields a ~ 0.5. For the reasons given above, 
we expect that in EX Hya a larger fraction is of ellipsoidal 
origin and adopt a = 0.6 ... 1.0 as a plausible range. 

We determine the flux / of the secondary star from A e \\ 
and an expression for the relative modulation A e \\/ f '. To 
this end, we again refer to the numerical study of Froning 
et al. (1999) who carefully determined the modulation 
using the if-band gravity and limb darkening coefficients 



of a Roche lobe filling secondary star of T c s — 3000 K, 
similar to that in EX Hya. Their result can be scaled as 



AelW/H = 0.13 



sin 2 i 



(2) 



where /h is the equivalent spectral flux of a spherical star, 
q = M2/M1 is the mass ratio, i is the orbital inclination, 
the functional form of Eq.|^ is that of the dominant term 
in Kopal's analytic expression (Binnendijk 1974), and the 
index H refers to the ff-band. With i = 79° (Hellier et al. 
2000), q ~ 0.19 ± 0.04 (Hellier et al. 1987, Vande Putte et 
al. 2001), and a = 0.8 ± 0.2 we then have 



/h = 7.7 aA 49i H 



1 



q 



■ 1 ■ 

sin 1 



= (7.6 ± 2.0) A 49 , H 



(3) 



where the error accounts for the uncertainties in a, ^49^, 
and q. With Ajg.H from Tab.|[ we obtain the _ff-band 
flux of the secondary star in EX Hya as /h = (9.2 ± 2.4) x 
10~ 16 ergcm _2 s _1 A _1 or 8.0±2.1mJy, where most of the 
error is systematic in nature and represents the full ac- 
ceptable range, not a standard deviation. Using the nu- 
merical factor of Eq. ^ for all infrared bands, we obtain 
the fluxes plotted in Fig. f| (open circles, statistical errors 
only). Because of the wavelength dependence of limb dark- 
ening (e.g. Al-Naimiy 1978), the true spectral energy dis- 
tribution may be slightly redder (with a pivot point at 
H). This is consistent with the spectral distr ibution of the 
secondary star determined below (Sect. O). Of the total 
orbital mean if-band flux of 13.3 mJy, this interpretation 
then assigns 7.3 mJy to the secondary star. Its infrared 
magnitude is K = 12.431°' 39- We continue the discussion 
of the spectrum of the secondary star in Sect. 4.3 . 

Of the remaining mean X-band flux, we attribute 
1.2 mJy to the the accretion funnel and 4.8 mJy to the 
accretion disk and white dwarf (Sects. 4.2 and [4.4| ). Since 
about 1/2 of the funnel emission and about 1/5 of the 
disk are eclipsed, we expect a central depth of the partial 
eclipse of about 1.7mJy or 1.1 x 10 -16 ergcm _2 s _1 A _1 , 
equal to the observed depth within the uncertainties 



(Sect. 3.2) 



4. Spectral analysis 

4.1. Overall spectrum 

Figure [| (upper left panel) shows the mean overall spectral 
energy distribution of EX Hya from the ultraviolet to the 
infrared. The mean IUE and HST spectra agree in abso- 
lute flux (Table |2|) and, in the region of overlap, we show 
the less noisy HST spectrum. We have also schematically 
included the ORFEUS11 continuum of Mauche (1999) 
which agrees reasonably well with the HUT spectrum of 
Greeley et al. (1997). In judging the overall spectrum, 
one should note that the individual sections were not ob- 
served contemporaneously. In spite of this, an astound- 
ingly coherent view emerges. The spectrum clearly shows 
the Balmer, the Paschen, and the Brackett jumps in emis- 
sion and, hence, contains substantial optically thin com- 
ponents, as suspected already by Berriman et al. (1985) 
on the basis of the infrared colours of EX Hya. 
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part of the spectrum (dot- 
ted). 



4.2. Spectrum of the spin-modulated component 

The spin-modulated component has been extracted at all 
wavelengths and the amplitudes Aq-j of the broad bands 
and the line fluxes are listed in Tables || and ||. The dif- 
ference between the spectral fluxes at spin maximum and 
spin minimum, 2Aq-j{\), is depicted in the lower panel of 
Fig.||. In the long-wavelength IUE range, the modulation 
amplitude could not be measured because of the lack of 
time resolution. We have, therefore, represented this spec- 
tral section schematically by the dotted line which has the 
mean slope of the long-wavelength spectrum and is ad- 
justed to the modulated flux of the short wavelength IUE 
spectrum. The modulation of the ORFEUSn spectrum 
at 1010 A (Mauche 1999) is indicated by a single data 
point. The spin-modulated component displays emission 
lines of the Balmer and Paschen series and the associ- 
ated jumps in emission, while the Brackett lines appear 
in absorption (see above). The mean full base width of 
the ultraviolet emission lines corresponds to velocities of 
±2300 km s" 1 (see Fig.|[ below), typical of the funnel at a 
couple of white dwarf radii. 

Spin modulation originates from the rotating partially 
optically thick accretion funnels at both poles of the white 
dwarf. Because any optically thin and non-occulted com- 



ponent is constant in time, the mean funnel flux /f un (A) > 
Aer(X). We need, therefore, an additional assumption to 
deduce /f un (A) from Aq7 (A). We note that the bright ul- 
traviolet lines of Cm, Civ, Sim, Siiv, and Nv are strongly 
modulated (see in Tab.||). For definiteness, we refer 
to the 80% modulation of the strongest ultraviolet line, 
CivA1550, and assume that this line originates entirely in 
the funnel which is not unreasonable because, as shown be- 
low, Civ absorption is weak in the white dwarf spectrum. 
Adopting the modulation of CivA1550 as representative of 
/fun (A), we compute the mean funnel flux to a first-order 
approximation as 



/ fun (A) = 1.5A e7 (X). 



(4) 



In the optical continuum, the eclipsed flux at spin maxi- 
mum slightly exceeds the spin modulated flux which sug- 
gests a similar factor. Nevertheless, there are clearly un- 
certainties in defining /f un (A) by a simple relation as in 
Eq. ||. We expect the ultraviolet emission of the heated 
polar caps of the white dwarf to display a spin modula- 
tion, too, and separation of the two sp in-m odulated com- 
ponents becomes non- trivial (see Sect. [4.5| ). On the other 
hand, /f un (A) clearly becomes small at long wavelengths 
and the implied error in the unmodulated component be- 
comes then small, too. We use Eq. ^ to compute the con- 



8 



S. Eisenbart et al.: Intermediate Polar EX Hydrae 





M3.0...M5.5 



5000 



10000 15000 20000 

Wavelength (A) 



Fig. 4. Contribution of 
the secondary to the mean 
spectrum of EXHya. The 
sample M-star spectra of 
spectral types M3.0, M4.5 
and M5.5, respectively, have 
been scaled to the expected 
fluxes (crosses) derived from 
the ellipsoidal modulations 
listed in Table || (see text). 
For clarity, the noisy region 
between the J and iJ-bands 
has been smoothed. 



tribution of the funnel emission to the luminosity, but do 
not attempt to model its spectral distribution. 

4.3. Spectrum of the secondary star 

Fig.|] shows the quasi- Johnson IJHK fluxes of the sec- 
ondary star as defined above (Sect. 3.2.2.) along with the 
spectra of M3 to M5.5 stars least squares fitted to the data 
points for a = 0.8. The systematic uncertainty in the flux 
is at most 2mJy (Sect. p~3| ). 

Formally, the four data points are best fitted by an 
M3 star, but, as noted above, the true spectrum will be 
slightly redder and the spectral type later, about M3 to 
M5, still consistent with the M3 assignment by Dhillon 
et al. (1997). The contribution of the secondary star to 
the total spectrum is largest at long wavelengths and the 
only pronounced distinctive feature in this range is the 
broad hump which is centered at ~ 2.2 fim and flanked 
by H2O absorption dips. This feature is clearly visible 
also in the spectrum of EX Hya. We determine the most 
probable contribution from the secondary star by sub- 
tracting the spectrum of an M3, an M4.5, and an M5.5 
star from the observed spectrum. Fig.|s| shows the corre- 
sponding difference spectra. We expect this spectrum to 
decrease smoothly since the main contributions should be 
from the spin-modulated component (Fig.^|) and the ac- 
cretion disk (Fig. [?]) and both decrease with wavelength. 
The smoothest decrease is obtained for the M4.5 star, in 
agreement with the expected spectral type in a CV with 
98 min orbital period (Beuermann et al. 1998). We adopt 
a spectral type M4 ± 1. 

The detection of absorption lines from the secondary 
in EX Hya is difficult. Fig.^J compares the spectrum of 
the M4.5 dwarf G1285 with spectra of EX Hya collected 
near spin minimum: (i) mean spin minimum, (ii) ellip- 
soidal maximum (</>gs ~ 0.75) near spin minimum, and (iii) 
eclipse (4>98 — 0, backside of star) near spin minimum. The 



M-star features are largely obscured by lthe emission line 
components in the 8300-9600 A range. In the if -band, the 
presence of the M-star is indicated by the faint Nal doublet 
and the CO molecular bands( see also Dhillon et al. 1997). 
Although our spectral resolution is not sufficient to resolve 
the Nal lines, phase-shifted superposition of the spectra 
suggests that Nal absorption is preferentially present on 
the dark side of the star and displays a radial velocity 
consistent with the motion of the secondary star (an un- 
certain 400 km s -1 ). Higher resolution is needed to settle 
the question of the radial velocity of the secondary star 
(Vande Putte et al. 2001). 

4.4. Spectrum of the accretion disk 

Figure shows the (slightly smoothed) spectrum which 
remains after subtraction of the funnel component /f un (A) 
and the secondary star, assuming it has K = 12.43, 
and is a dM4.5 star (Figs, [l] and ||). We assign this 
spectrum to the sum of the accretion disk including the 
bright spot at its edge and the white dwarf. The spec- 
trum lacks the intense ultraviolet lines which belong to 
the funnel component, but contains a substantial fraction 
of the Hei, Balmer, Paschen, and Brackett line emission. 
MgnA2800 may, in part, originate in the funnel and be 
falsely assigned to the disk because we subtracted the 
long-wavelength IUE contribution to the funnel spectrum 
only in a schematic way without accounting for the emis- 
sion lines. The (unsmoothed) short-wavelength section of 
the spectrum in Fig. |?] is shown enlarged in the insert and 
is further discussed in Sect. [D] (Fig. ||). This part of the 
spectrum is dominated by the heated white dwarf. We 
model the A > 3000 A part of the spectrum by an isobaric 
and isothermal disk of pure hydrogen plus a blackbody 
representation of the white dwarf. The spectrum of the 
disk is computed following the description of Gansicke et 
al. (1997). It has an outer radius of 1.6 x 10 10 cm, a thick- 
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Fig. 6. Spectral signatures of the secondary star in the IR spectrum of EXHya. Both panels show the spectrum 
of the M4.5 star G1285 and three spectra of EXHya, (i) mean near spin minimum, (ii) 49-min maximum near spin 
minimum, and (iii) eclipse near spin minimum (from top to bottom). All spectra are normalised to unity at 9800 A (left 
panel) or 22500 A (right panel) and vertically offset by multiples of 0.5 flux units (left panel) and 0.3 flux units (right 
panel). Spectral features expected to be present in the secondary star are indicated. H2O refers to regions affected 
by atmospheric water vapour absorption. The strong emission lines in the EX Hya spectra are the hydrogen Paschen 
series and CaiiA8498, 8542, 8662 (left panel) and B7 (right panel). 



ness of 2 x 10 8 cm, a pressure of 200 dyne cm~ 2 , an electron 
temperature of 7000 K, close to the canonical tempera- 
ture of an optically thin accretion disk (Williams 1980), 
and is seen at an inclination of 79°. For definiteness, we 
assume that the disk has a central hole of 6 x 10 9 cm ra- 



dius (Heflier et al. 1987), but the structure could also be 
a ring with a larger inner radius surrounding the magne- 
tosphere of the white dwarf. This component is optically 
thin and accounts for much of the line emission and the 
Balmer, Paschen, and Brackett jumps, but lacks contin- 
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uum flux around 5000 A. We include, therefore, a second 
component which has a pressure of 10 3 dyne cm~ 2 , an elec- 
tron temperature of 10 000 K, a weak Balmer jump, and is 
nearly optically thick. The hotter component covers only 
~ 20% of the area of the cool one and may r epre sent the 
heated inner edge of the disk (but see Sect. 4.5). Fig. || 
shows the sum of the two disk components (short-dashed 
curve). Although this is not a formal fit, it reproduces the 
observed optical and infrared part of the spectrum reason- 
ably well. Metal line emission (mostly Fen, not included in 
the model) may substantially contribute to the observed 
excess in the range 4000. . . 5500 A. Lines of the Pfund se- 
ries (not included either) may contribute to the excess flux 
longward of 22 788 A. 

The spectrum in Fig. [?] is based on K = 12.43 for 
the secondary star. If the secondary were substantially 
brighter, the decreased infrared flux could not be fit by 
the optically thin spectrum anymore; if it were fainter, 
there could be additional disk flux ove r that depicted in 
Fig. within the limits set in Sect. The 2 mJy uncer- 
tainty in the iVband flux of the secondary, would allow 
an underlying cool component of the disk with a black- 
body temperature up to Tbb = 3300 K for a disk radius 
of 1.6 x 10 10 cm. This is a significant result because the 
mass in the modeled disk is <, 6 x 10 19 g, much less than 
the ~ 10 22 g needed to feed one of the rare outbursts. If 
the outbursts result from a disk instability, a layered disk 
is needed in which the component of high surface density 
must be cooler than the limit given above. This result is 
very similar to the conclusions of Gansicke et al. (1997) 
for the case of the SU UMa dwarf nova EK Tra. 

Fig. ^| shows the HST spectrum before and after 
eclipse, at P rb = ±(0.05 — 0.15). Before eclipse pro- 



nounced selfabsorption cores of highly excited ultravio- 
let lines are seen which disappear after the eclipse. These 
cores are produced in the bulge at the edge of the ac- 
cretion disk, i.e. the interaction region of the accretion 
stream with the disk or ring. This bulge produces also the 
enhanced Paschen line emission at (f>gg, = 0.25, the orbital 
maximum in the i?-band at </>g§ = 0.85 (Hellier et al. 2000), 
and the pre-eclipse dip in the EUV and soft X-ray regime 
(Cordova et al. 1985, Rosen et al. 1988, Hurwitz et al. 
1997, Mauche 1999). The radial velocity of the absorption 
lines in the difference spectrum (C in Fig. ||) is consistent 
with zero, as expected for material at the edge of the disk 
or ring which moves essentially perpendicular to the line 
of sight. The presence of Nv indicates an excitation tem- 
perature of the order of 50 000 K or higher. The metal lines 
are just resolved with a FWHM of 2.5 A, corresponding to 
about 400 km s . The Lya absorption line in the differ- 
ence spectrum is wider with FWHM ~ 5 A (1000 kms" 1 ), 
suggesting an origin in warm material or material with a 
large velocity dispersion. The column density needed to 
produce the bulge absorption in Lya at (/>gs~ 0.9 is about 
10 19 H-atomscm -2 , substantially less than the equivalent 
column density of 1.3 x 10 20 H-atomscm -2 of cold mat- 
ter of solar composition needed to account for the ob- 
served EUV bulge absorption (Hurwitz et al. 1997). The 
difference suggests that hydrogen is largely ionized and 
the EUV bulge absorption is primarily due to helium. 

4.5. Spectrum of the white dwarf 

The ultraviolet section of the spectrum in Fig.[?] 
is reproduced in Fig. |^. It shows the typical ab- 
sorption lines expected from a hot atmosphere, in- 
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eluding CiiiA1175, SiiiA1260, 1265, Sin, mA1294 ... 1306, 
CiiA1335, SiivA1394, 1403, and SinA1527, 1533. Some 
of the lines are affected by the subtraction proce- 
dure by which this spectrum was created, CivA1550, 
SiivA1394, 1403, and probably ClllA1175. In order to es- 
timate the temperature of the white dwarf (or its heated 
pole cap) in EXHya, we have computed a set of white 
dwarf spectra covering T c g — 15 000 — 50 000 K in steps of 
5000 K, using the codes TLUSTY 195 and SYNSPEC 45 
(Hubeny 1988; Hubeny & Lanz 1995), assuming solar 
abundances and log g = 8.0 (i.e. M Q ~ 0.6). The ob- 
served line ratios Sim/Sin and Cm/Cii constrain the 
temperature of the white dwarf photosphere to T G ff = 
25 000 ± 3000 K and a white dwarf model of this temper- 
ature qualitatively fits also the absorption lines of other 
heavy elements. For an distance of 65 pc (Sect. |5.2j ), the 
observed flux in the range 1300 - 1500 A corresponds to a 
white dwarf radius of 6 x 10 s cm. 

We do not attempt a detailed fit to the white dwarf 
spectrum for the following reasons: (i) the irradiated white 
dwarf has very likely a non-uniform temperature distribu- 
tion; (ii) the spectrum is probably altered by attenuation 
in the accretion curtain; and (iii) the long wavelength part 
of the white dwarf spectrum is heavily veiled by the ac- 
cretion disk and practically unobservable. We comment on 
these points in turn. 

The funnel radiation which impinges on the white 
dwarf will produce extended hot polar caps. Taking ac- 
count of these effects would introduce additional parame- 
ters which are not well constrained. Furthermore, heating 
affects the temperature structure of the white dwarf at- 
mosphere and thereby alters the profiles of the absorption 
lines. In the polar AM Her, where the FUV emission is 
dominated by emission from the white dwarf, we could 
show that irradiation results in a significantly decreased 
depth of the Lya absorption line (Gansicke et al. 1998). 



A similar effect could explain the shallow Lya absorption 
observed in EX HYa. 

The mean photospheric spectrum of the white dwarf is 
probably altered by the fact that part of the light is atten- 
uated on its passage through the accretion curtain. An in- 
dication for this is provided by the double hump structure 
of the spectrum shown in the insert in Fig. |?]. Such spectral 
shape is reminiscent of the white dwarf spectrum in OY 
Car which is heavily veiled by metal absorption in material 
referred to as "Fell curtain" (Home et al. 1994). There is, 
in fact, a detailed similarity between the ultraviolet spec- 
trum in Fig. and that of the white dwarf in OY Car (Fig. 
5, bottom panel of Home et al.). Metal absorption in the 
funnel produces the broad dips around 1700 A and 2500 A, 
and the remnant spike at 1595 A. Given that funnel emis- 
sion and absorption is abundant at other wavelengths, we 
are not surprised to find funnel absorption features mixed 
up with photospheric ones. The strong absorption lines in- 
dicated in Fig. ^ are not noticeably broadened, however, 
and probably of photospheric origin. 

If part of the ultraviolet emission is due to a white 
dwarf with heated polar caps, we may expect that the flux 
shows spin-modulation. There are spin-modulated contri- 
butions to the ultraviolet flux by the accretion funnel and 
by the heated white dwarf. The best visibility of the ac- 
cretion funnel is when the lower pole points towards the 
observer (Beuermann & Osborne 1988). At this rotation 
phase, </>67 = 0, the heated polar caps reach their mini- 
mum projected area. Hence, the spin modulated compo- 
nent from the heated polar caps assumes its minimum 
when the funnel flux reaches its maximum, i.e., the two 
modulations will be out of phase by 180°. The amplitude 
of the spin modulation will then display a minimum near 
1500 A, where the spectral flux of the white dwarf reaches 
a maximum, as it is observed (Fig.|^). 
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Finally, there will be a contribution from the unheated 
fraction of the white dwarf photosphere. For an effective 
temperature of 10 000-15 000 K and 10 9 cm radius, the 
peak flux would reach only 1-3 mJy and can not be identi- 
fied against the other components. It might be contained, 
e.g., in the adopted disk flux which would then have to be 
reduced correspondingly (Sect. p^ ). 

5. System parameters 

5.1. Mass and radius of the secondary star 

We assume that the secondary is a main sequence star 
which is negligibly expanded over its zero age main se- 
quence radius (Baraffe et al. 1998, Kolb & Baraffe 1999, 
see also Beuermann et al. 1998). For stars of solar compo- 
sition, the Baraffe et al. mass radius relation of spherical 
stars may be approximated by a power law of the type 
R oc M@ for limited ranges of mass. Considering that the 
equilibrium radius of a Roche lobe filling star of the same 
mass is larger by about 6% (Kolb 2001), we obtain 

R 2 /R Q =0.77 (M 2 /M ) ' 75 (5) 

for masses between 0.11 and 0.40 Mq. The relation of 
Patterson (1984), R 2 /R G = (M 2 /M©) - 88 , predicts radii 
within 5% of this for masses between 0.09 and 0.20 Mq. 
Eq. (§) and Roche geometry yield M 2 = 0.12 ± 0.01 M© 
and R 2 = (0.157 ± 0.010) R Q = (1.09 ± 0.07) x 10 10 cm, 
where we have allowed for some systematic error in Mi. 
If the mass exceeds 0.13 M©, the secondary would have 
to be smaller than a main sequence star. If it is as low as 
0.10 Mq, an expansion in radius by more than 10% would 
be required. 

5.2. Distance of EX Hydrae 

The detection of the secondary star in EX Hya allows 
to obtain an estimate of the distance d by the surface 
brightness method (Bailey 1981), 

\ogd= (K-S K )/5 + l + \og(R/R Q ). (6) 

We use the re-calibration of the i-T-band surface brightness 
of late type stars by Beuermann (2000) which yields 5k = 
4.37±0.22 for a dwarf of spectral type M3-5 and near-solar 
composition. With K = 12.19 . . . 12.73 and the stellar ra- 
dius of the last section, we obtain a distance d = 65±llpc. 
This is much lower than previous estimates which include 
105 pc (Warner 1987) and > 130 pc (Berriman et al. 1985). 
EX Hya may be one of the closest CVs with a known dis- 
tance. 

The largest contribution to the error in d arises from 
the uncertainty in the ^-magnitude of the secondary 
(Sect. pT3| ). The distance is near the upper limit of the 
derived range if the true ellipsoidal modulation is as small 
as 60% of the observed 49-min variation. Because of this 
uncertainty and the special importance of EX Hya among 
all CVs, a trigonometric parallax for this object is needed. 



Table 4. Energy fluxes of the spectral components of 
EX Hya. 



Component 


hi (kevj 


Flux 




or A (A) (10" 


11 erg cm _2 s _1 ) 


Hard X-rays 


> 1.0 keV 


15 


Soft X-rays 


0.28...1.0keV 


3 


EUV 


0.067.. . 0.28 keV 


~ 1 




0.0136.. . 0.067 keV 


~ 1 : 


Accretion curtain 


A = 912. ..24 000 A 


12 


Accretion disk 


A = 912. ..24 000 A 


10 


White dwarf 


A = 912. ..24 000 A 


18 


Total 




60 



5.3. Interstellar absorption 

X-ray spectral fits suggest that the X-ray source is par- 
tially covered by substantial column densities TVh of hy- 
drogen (e.g. Allan et al. 1998). This is material within the 
binary system. The detection of EUV emission (Hurwitz 
et al. 1997) proves that the true interstellar column den- 
sity is small. Fitting the narrow Lya-absorption profile 
in the 0gg = 0.05-0.15 spectrum of Fig. || yields Vh — 
(3 ± 1) x 10 18 H-atoms cm -2 , consistent with the non- 
detectable extinction, Eb-v 0.05 (Verbunt 1987). At 
d = 65 pc, the mean space density of neutral atomic hydro- 
gen along the line of sight to EX Hya (I, b = 304°, +33°) 
is n H = 0.01... 0.02 cm" 3 . 

5.4. Luminosity 

We estimate the luminosity as L = 4 ir d 2 F, where F is the 
orbital and spin-averaged energy flux, with the secondary 
star subtracted and corrected for interstellar absorption 
using Ah = 3 x 10 18 H-atoms cm -2 of cold matter. There 
is some uncertainty in the estimate of L because the mean 
flux used may deviate from the 47r-average. 

Table ^] lists the mean fluxes in individual energy 
bands shortward of Lya and for individual physical com- 
ponents longward of Lya. The mean X-ray flux is col- 
lected from Beuermann & Osborne (1984), Cordova et al. 
(1985), Rosen et al. (1988), Rosen et al. (1991), Singh 
& Swank (1993) and Allan et al. (1998). The observed 
EUV and soft X-ray flux in the 0.067-0.280 keV range is 
~ 0.7 x 10~ n erg cm~ 2 s _1 (Cordova et al. 1985, Hurwitz et 
al. 1997) which after correction for interstellar absorption 
increases to ~ 10 erg cm s . There is no informa- 
tion on the flux between 13.6 and 67 eV and we simply 
add another ~ 10 -11 ergcm _2 s _1 . 

Angle-dependent absorption and scattering in the fun- 
nel is probably responsible for the observed low-energy 
depression and the spin modulation of the EUV and X- 
ray spectrum (e.g. Kim & Beuermann 1995). The flux 
contained in the spin-modulated soft and hard X-ray 
components, i.e. half the difference between maximum 
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and minimum fluxes for photon energies > 0.28 keV, 
amounts to 3 x 1CP 11 ergcm~ 2 s -1 . Applying the Zanstra 
method to the flux in the HeiiA4686 line, i^se = 
1.65 x 10~ 13 ergcm _2 s _1 (Patterson & Raymond 1985, 
this work), demonstrates that an additional flux -F>54 ~ 
3 x 10~ n ergcm -2 s -1 is absorbed at photon energies be- 
tween 0.054 and 0.28 keV. The sum of these two com- 
ponents is the minimum reprocessed flux expected from 
the accretion funnel. The funnel emission will actually 
be higher because additional energy is absorbed from 
the Lyct continuum below 54 eV and from the Balmcr 
continuum of the heated white dwarf. Furthermore, vis- 
cous interaction of the accreted gas in the funnel may 
release additional energy. The sum of all these sources 
should account for the mean observed funnel emission of 
12 x 10 -11 ergcm _2 s _1 which appears reasonable. The lat- 
ter number is the wavelength-integral of the funnel com- 
ponent for A > 912 A , F fun = 1.5 J A 67 dX (see Sect. |t] 
and Fig. g). 

Subtracting the funnel component and the secondary 
star from the mean spectrum leaves the accretion disk and 
the white dwarf (Fig. ^). In Table [|, the disk is represented 
by the integral over the dotted spectrum in Fig. [?] and the 
white dwarf by the flux remaining after subtraction of the 
dotted spectrum. The disk component contains about 1/6 
of the total flux (Table |^) a bit more than expected from 
gravitational energy release outside r\ ~ 6 x 10 9 cm and a 
white dwarf of radius R\ — 10 9 cm. Note that this result 
does not allow an inner disk radius n < 6 R\ , consistent 
with the visibility of the lower pole of the white dwarf in 
X-rays (Beuermann & Osborne 1988). 

The intrinsic flux of the white dwarf can not be sepa- 
rately identified and is included in the contribution quoted 
for either the heated white dwarf or the accretion disk 
(Table ||). For R\ — 10 9 cm and an effective tempera- 
ture of the white dwarf of 10 000-15 000 K, it amounts 
to only (2 — 6) x 10~ n ergcm^ 2 s _1 . Hence, the total 
accretion-induced flux is slighly reduced below the sum 
of the fluxes in Table [| to about 5.6 x 10 -10 ergcm~ 2 s -1 . 
This corresponds to a luminosity of L = 2.8 x 10 32 ergs -1 
at d = 65 pc, or L = 6.7 x 10 32 d 2 00 ergs -1 , where 
g?ioo is the distance in units of 100 pc. Not included in 
Table ^ is the secondary star which contributes only 
1.5 x 10 -11 ergcm~ 2 s -1 . 

5.5. Accretion rate 

The accretion rate is obtained as M = LR±/ GM\ , 
where M\ and R\ are mass and radius of the white 
dwarf. With L from the last section, M = 5 x 10 15 
( J Ri/10 9 cm)(A/i/A/ Q ) -1 (i 1 2 00 gs -:L . M x is estimated from 
the X-ray temperature as 0.45 - 0.48 M Q (Fugimoto & 
Ishida 1997, Cropper et al. 1999, Ezuka & Ishida 1999, 
Ramsay 2000). The radial velocity amplitudes of the white 
dwarf, K x ~ 69 ± 9kms -1 (Hellier et al. 1987), and of 
the secondary star, Kz — 356 ± 44kms _1 (Vande Putte 
et al. 2001), yield M 2 = 0.095 M and Mi = 0.49 Mq. 



According to our discussion in Sect. |5.l[ the secondary 
would have to be substantially expanded over a main se- 
quence star if of such a low mass and we prefer the higher 
mass of 0.12 M Q because, within the errors, K\ and K2 
are consistent also with M2 = 0.12 Mq and Mi as high as 
0.7 Mq. Given the remaining uncertainty in the masses, 
we quote the accretion rate in terms of a standard white 
dwarf of 0.6 Mq. Using a Wood (1994) white dwarf model 
with thick hydrogen envelope and T c g— 10 4 K, the radius 
is then Ri = 8.9 x 10 s cm and M ~ 7.5 x 10 15 d? m gs -1 or 
M ~ 3.1 x 10 15 gs -1 at d = 65 pc. This is consistent with 
the accretion rate as expected from gravitational radiation 
as the dominant angular momentum loss mechanism. It 
is also consistent with the distance-independent accretion 
rate deduced from the observed spin-up rate of the white 
dwarf in EX Hya, M spin = 2.4x 10 15 (n/10 10 cm)" 1 / 2 gs" 1 
for Mi = 0.6 Mq (Ritter 1985) which is derived on the as- 
sumption that the accreted matter couples onto the mag- 
netic field of the white dwarf at the inner disk radius r; . 
The estimated n ~ 6 x 10 9 cm (Hellier et al. 1987) yields 
M sp i n — 3.1 x 10 15 gs -1 . These numbers do not change 
substantially if M\ ~ 0.5 Mq. 

5. 6. Magnetic field strength of the white dwarf 

A coarse estimate of the (di-)polar field strength of the 
white dwarf of B p ~ 0.18 MG can be derived from n ~ 
6 x 10 9 cm (Hellier et al. 1987) and M~3x lO^gs -1 if 
n is equated to r M ~ 2~ 10 / 7 fi 4 / 7 (GM^- 1 / 7 M~ 2 / 7 cm, i.e. 
about 1/2 of the spherical Alfven radius (e.g. Frank et al. 
1992, their Eq.6.12). Larger r; (Hellier et al. 2000) would 
lower the field estimate. Such a low field is consistent with 
the lack of cyclotron emission in the infrared. 

6. Conclusions 

We have detected the dM4 ± 1 secondary star in EX Hya 
by its ellipsoidal modulation and spectral energy distri- 
bution. From its infrared magnitude K ~ 12.5, we have 
derived a distance of EX Hya of 65 ± 11 pc, substantially 
less than previous estimates. The integrated flux from the 
X-ray to the infrared range averaged over spin and orbital 
period is 6 x 10~ 10 ergcm -2 s -1 , yielding a luminosity of 
about 3 x 10 32 ergs -1 and an accretion rate which is con- 
sistent with the observed spin-up and with gravitational 
radiation as the dominant angular momentum loss mech- 
anism. Distance, luminosity, and accretion rate depend 
upon which fraction of the observed infrared continuum 
modulation at half the orbital period may be interpreted 
as ellipsoidal modulation of the secondary star. If this frac- 
tion is near unity, d is near the lower limit of 55 pc, if it 
less than 60%, d may exceed 75 pc. This question can best 
be resolved by obtaining a trigonometric parallax. 

The heated white dwarf is detected for the first time by 
strong metal absorption lines which are probably of pho- 
tospheric origin, but may be affected also by absorption 
in the accretion funnel. The equivalent width ratios of the 
Sin, in, iv and C11, ill lines suggest a near-solar metal 
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abundance in the uppermost photosphcric layers and an 
effective temperature near 25 000 K. The radius needed 
to account for the emission is about 6 x 10 8 cm. If the 
white dwarf is of low mass as suggested by several stud- 
ies and has a radius of about 10 9 cm then what we see 
are probably the heated polar caps of the white dwarf. 
Intrinsically, the white dwarf may be substantially cooler 
than indicated by the line ratios (and perhaps display a 
less metal- rich surface composition) . The observed "white 
dwarf flux" corresponds to ~ 28% of the total accretion 
flux, consistent with the reprocessed fraction of the down- 
ward flux from the post-shock region considering both, 
geometry and reflection probability. 

There is agreement that maximum light at optical 
and X-ray wavelengths occurs when the upper pole points 
away from the observer (Hellier et al. 1987, Beuermann & 
Osborne 1988, Rosen et al. 1988, Siegel et al. 1989, Hellier 
et al. 2000). At this spin phase, the light from the heated 
caps is likely to reach a minimum which implies antiphased 
modulations of the ultraviolet light from the funnel and 
the white dwarf surface. The minimum in the amplitude 
of the spin modulation near 1500 A is consistent with this 
notation. More extensive ultraviolet observations can pos- 
sibly disentangle the contributions from both components. 

Of the outward directed integrated X-ray flux, the 
spin-modulated fraction of ~ 20% is absorbed and re- 
processed in the accretion funnel. The funnel spectrum is 
largely the result of radiative transfer of the X-ray emis- 
sion from the accretion column and of the ultraviolet emis- 
sion from the heated white dwarf. Although it is clear 
that the funnel is optically thick in the emission lines and 
partly in the continuum, the detailed funnel geometry re- 
mains uncertain and we have not attempted to model the 
resulting spectrum (Fig. 0). Better spectral resolution in 
the optical and infrared is needed to study the radiative 
transfer in detail. 

The emission which is left after subtraction of the sec- 
ondary star, the white dwarf, and the funnel emission 
originates from the accretion disk or ring surrounding the 
white dwarf magnetosphere. The spectral flux contribu- 
tion of much of the remaining optical and infrared light is 
plausibly explained by an isothermal slab of gas of radius 
1.6 x 10 10 cm, thickness 2 x 10 8 cm, and a temperature of 
most of the gas near 7000 K which is optically thin plus a 
largely optically thick component at 10 000 K which cov- 
ers only 20% of the emitting area of the disk. There is 
a tight upper limit on the blackbody temperature of an 
underlying optically thick component covering the entire 
disk of T bb < 3300 K. 

In summary, we have presented an interpretation of 
the overall spectral energy distribution of EX Hya which 
explains the individual components within the standard 
model of a weakly magnetic intermediate polar. 
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